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In this paper we explore the connection between anisotropic Gaussian fluctuations and isotropic 
non-Gaussian fluctuations. We first set up a large angle framework for characterizing non-Gaussian 
fluctuations: large angle non-Gaussian spectra. We then consider anisotropic Gaussian fluctuations 
in two different situations. Firstly we look at anisotropic space-times and propose a prescription 
for superimposed Gaussian fluctuations; we argue against accidental symmetry in the fluctuations 
and that therefore the fluctuations should be anisotropic. We show how these fluctuations display 
previously known non-Gaussian effects both in the angular power spectrum and in non-Gaussian 
spectra. Secondly we consider the anisotropic Grischuk-Zel'dovich effect. We construct a flat space 
£ — ■ time with anisotropic, non-trivial topology and show how Gaussian fluctuations in such a space- 

^\ | time look non-Gaussian. In particular we show how non-Gaussian spectra may probe superhorizon 

■ anisotropy. 
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PACS Numbers : 98.80.Cq, 98.70.Vc, 98.80. Hw 

I. INTRODUCTION 



Anisotropic models of the Universe have been often considered in the past (eg. [|l]]2|]). In recent times globally 
anisotropic spacetimes have attracted attention for their thought provoking value, as primordial anisotropy would 
appear to contradict inflation It is therefore important to find experimental evidence for, or constraints on, 
primordial anisotropy. The cosmic microwave background (CMB) is the cleanest and most accurate experimental 
probe in current cosmology. Thus it makes sense to explore the impact of anisotropic expansion on the CMB. For 
homogeneous space times this was largely done in [^-^]. In the more sophisticated analysis in || the effects of the 
unperturbed anisotropic expansion were combined with a spectrum of superposed Gaussian fluctuations. An admitted 
shortcoming of this analysis is the assumption that while the unperturbed model leaves an anisotropic pattern in the 
' sky, the Gaussian fluctuations around it are isotropic. Should the Gaussian fluctuations in such models be anisotropic 
one may expect a more stringent statistical bound on anisotropy, if the Universe is indeed isotropic. One can consider 
O another class of models where the background space-time is homogeneous and isotropic but anisotropic topological 
identifications lead to anisotropic Gaussian fluctuations. Some of these universes have been considered before |?J and 
an example of the patterns in an open universe has been presented in ||] . 

The apparently unrelated issue of large-angle CMB non-Gaussianity has also been considered recently, both as an 
experimental matter ||, and as a possible prediction in topological defect theories |lfj|-|Tq|. In jHJ, in particular, 
an outline is given of a comprehensive formalism for encoding large angle non-Gaussianity based on the spherical 
• i-H , harmonic coefficients a l m in the expansion 

t=0m=-t 

In |To| it is also stated that "any non-Gaussian theory is to some extent anisotropic, favouring particular directions 
in the sky and some m's over others" . The converse statement follows: that Gaussian anisotropic fluctuations will 
appear as non-Gaussian fluctuations from the standpoint of an isotropic theory. This establishes an interesting link 
between the search for cosmological anisotropy and the search for non-Gaussian signatures. 

Let us consider Gaussian theories which favour an axis, where fi are angles defining this axis. Then the probability 
distribution conditional to this axis P(al n \il) is Gaussian. Isotropy is violated, but the resulting theory is Gaussian 
within the reduced set of symmetries the theory now must satisfy. However from an isotropic point of view the full 
ensemble is made up of all the ensembles which favour an axis, but allowing the axis to be uniformly distributed. Such 
a super-ensemble would undoubtedly be isotropic, but it would also be non-Gaussian. Marginalizing with respect to 
the axis reveals a non-Gaussian theory, that is 



P(<4) = J dQP(Sl)P(ai\n) 



P(O) = =! (2) 
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is non-Gaussian. This identifies the origin of the Gaussian/non-Gaussian switch. Conditionalizing to an axis renders 
the theory Gaussian (and anisotropic). Marginalizing with respect to the axis reveals a non-Gaussian theory (but an 
isotropic ensemble). 

This phenomenon turns out to be a particular case of the general phenomenon discussed in connection with the 
texture analytical model in pl|JT^ ]. In that model it is found that the temperature anisotropies are very non-Gaussian. 
The theory has Ci cosmic variance error bars above their Gaussian value, and there are strong correlations among 
Ce. It turns out, however, that these large-angle non-Gaussian effects are largely due to the last texture (as in the 
texture closest to us, or the texture at lower redshift). The culprit identified, one then notices that conditionalizing 
the theory to the last texture redshift z\ reveals a Gaussian ensemble, that is, the probability distribution P(a l m \z\) 
is Gaussian. Marginalizing with respect to z\, however, produces a non-Gaussian ensemble, that is the probability 

P{a l m ) = J dzxP^Plptizi) (3) 

is non-Gaussian. 

The picture is then clear fl7]| . We come up with a construction where the full ensemble is made up of sub-ensembles 
which are Gaussian. Each sub-ensemble is however labelled by an index which from the the point of view of the 
full ensemble is a random variable. Marginalizing with respect to this variable reveals a non-Gaussian ensemble. 
Conditionalizing with respect to this index renders the theory Gaussian. Such an index was called in |Q | the random 
index, and it was conjectured [] in that paper that non-Gaussianity could often be characterized by a set of such indices 
labelling Gaussian ensembles. Within such a construction the strategy for predicting experiment must be modified. 
One should now not provide a direct statistical description of the full ensemble (that is, marginal distributions), which 
would be plagued by all sorts of non-Gaussian effects. Rather it makes more sense to supply information on all the 
Gaussian sub-ensembles, plus the distribution function of their random indices. 

Hence we may use a sub-class of the comprehensive formalism for encoding large-angle non-Gaussianity outlined 
in |h]] to describe anisotropic Gaussian fluctuations. This is essentially a large-angle generalization of and is 
described in Section ||. The idea is to complement the angular power spectrum Ci with a set of multipole shape 
spectra B m describing how the power is distributed among the m's for a given scale I. The B m encode information 
on the shape of large angle structures. They are uniformly distributed in a Gaussian isotropic theory, meaning its 
fluctuations are shapeless. However, as we shall see in Section preferred shapes emerge in non-Gaussian isotropic 
theories, as well as in Gaussian anisotropic theories, where the B m are not uniformly distributed. Non-Gaussian 
spectra then appear as a natural predictive tool for these theories. 

In this paper we study the disguised non-Gaussianity of anisotropic Gaussian fluctuations along two lines. Firstly, in 



Section III, we propose a simple method for defining anisotropic Gaussian fluctuations. Breaking isotropy essentially 
amounts to choosing an alternative symmetry group under which the covariance matrix should be invariant, and which 
picks a favoured direction in the sky. We can then write down the most general form for the covariance matrix of the 
theory simply by studying the representation theory of the symmetry group. We argue that the accidental symmetry 
allowing anisotropic fluctuations to be isotropic is a model dependent and unnatural assumption. Hence Gaussian 
fluctuations in anisotropic Universes should be anisotropic too. Although we concentrate on anisotropic fluctuations 
with an SO(2) symmetry, the definition and considerations given in Section II arc quite general, as explained in more 
detail in Appendix I. 

We then show how anisotropic Gaussian theories induce well known non-Gaussian effects in the relation between 
the observed and the predicted angular power spectrum Ce- These effects include larger cosmic variance error bars, 
and also the phenomenon of cosmic covariance, that is correlations between the observed Ci. Cosmic covariance 
allows for more structure to exist in each realization than in the predicted average power spectrum and complicates 
comparison between theory and experiment. These effects are shown to be present for anisotropic Gaussian theories 



in Section IV. 

Then, in Section [v| we show how anisotropic Gaussian fluctuations render non-Gaussian spectra non-uniformly 
distributed, as announced above. We also find the most general class of isotropic non-Gaussian theories into which 



anisotropic Gaussian fluctuations may be mapped. As a concrete example in Section VI we proceed to characterize 
the non-Gaussian spectra for the relevant , glo bally anisotropic space times. 



Along a totally different line in Section VII we construct a simple example of a topologically non-trivial space time 



and show how the non-Gaussian spectra will indicate anisotropic topological identifications. We propose this as an 



1 This conjecture can in fact be promoted to a mathematical theorem; see fljf . 
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anisotropic Grischuk-ZePDovich effect: from subhorizon, large angle observables we can characterize super-horizon 
anisotropies. 

In Section VIII we discuss the implications of our results and their practical implementation. 



II. LARGE-ANGLE NON-GAUSSIANITY 



We now set up a formalism for describing large-angle non-Gaussianity which is based on [|19f| , but makes use of 
coefficients rather than Fourier components, and so is suitable for mapping large-angle non-Gaussianity. Again the 
idea is to map the {a^} into a set of spectra which for a Gaussian isotropic theory are independent random variables. 
One of these spectra is the angular power spectrum Cg, and should be a x\i+\ f° r a Gaussian isotropic theory. The 
other variables make up non-Gaussian spectra which should be uniformly distributed for a Gaussian isotropic theory. 

The transformation proposed is defined as follows. Firstly we split the complex modes into moduli and phases 

t 11 

a o — s oPo 

at = ^|e*« (4) 

where s l Q — ±1 is simply the sign of Oq. The fact that the m = mode is real introduces a slight modification to the 
construction in There are now £+1 moduli, but there are only £ phases (the index m starts at 1 for the phases). 
Working out the Jacobian of the transformation shows that for a Gaussian theory the distribution of the 4> m , Sg} 
is 



2exp 



\ 2Ct \ i i 

F(PL tin, 4) = ~ , 1/9 ^+l/2 >< II P™ P T^Al X 9 (5) 



(27T)V2^ +1/2 Vl J (27F) 



The phases (j) m are uniformly distributed in [0, 2tt\. The sign Sq has a uniform discrete distribution. The moduli 
p l m are x| distributed except for p l which is xl distributed. Since p now does not appear in the Jacobian of the 
transformation, the only way one can proceed with the construction in ]l9| is by ordering the p's by decreasing order 
of m, and then introduce polars: 

p\ = r cos 8i 
p e e _ 1 = rsin^i cos 02 

p{ = v sin cos 9i 

p = r sin#i...sin#; (6) 

Again, working out the jacobian of the transformation implies that for a Gaussian isotropic theory the distribution 
of these variables is 

F(r,e m ,sUm) - ^^^ncos^(sin^) 2 (-) x \ x JL (7 ) 
One can then define shape spectra as 

< = (sinft m ) 2 ^- m ) +1 (8) 
so that for a Gaussian isotropic theory one has: 

F(rB* sU*) ^(-r 2 /2C e yt 1 1 

Fir, B m , s , rfj - (7r/2)1/2c , +1/2(2£ _ 2 {27T y (9) 

The angular power spectrum C seen as a random variable is then related to r by 

r 2 

C l = (10) 

2£+l v ' 
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and is a Xm+v The multipole shape spectra may be obtained from the moduli p l m according to 

Bfn = ( P V+ (11) 

and are uniformly distributed in [0, 1]. Finally the phases are uniformly distributed in [0, 2w], and the sign Sg is 
a discrete uniform distribution over {— 1,+1}. 

As in Jl| we define non-Gaussian structure in terms of departures from uniformity and independence in the 
{B^jCj)^}. Gaussian theories can only allow for modulation, that is, a non-constant power spectrum. The most 
general power spectrum has as much information as Gaussian theories can carry. White-noise is the only type of 
fluctuations which is more limited in terms of structure than Gaussian fluctuations 0. In isotropic Gaussian theories 
there is no structure in the {B^, (j/ m } since these are independent and uniformly distributed. By allowing the B l m to 
be not uniformly distributed, or to be constrained by correlations amongst themselves and with the power spectrum, 
one adds shape to the multipolcs. This is because the B^ n tell us how the power in multipole I given by C e (or r) 
is distributed among the various \m\ modes, which reflect the shape of the fluctuations. Indeed the m = mode 
(zonal mode) has no azimuthal dependence. It corresponds to fluctuations with strict cylindrical symmetry (rather 
than statistical symmetry). The \m\ > modes correspond to the various azimuthal frequencies allowed for the scale 
I. Each of these modes represent a way in which strict cylindric symmetry may be broken. The relative intensities 
of all the to modes carry information on the shape of the random structures at least as seen by the scale I. In a 
Gaussian theory all the to modes must have the same intensity, something which can be rephrased by the statement 
that the B l m are independent and uniformly distributed. Hence Gaussian fluctuation display shapeless multipoles. 
Any departure from this distribution in the B^ may then be regarded as a evidence for more or less random shape 
in the fluctuations. 

On the other hand the phases (fr m transform under azimuthal rotations. Therefore they carry information on the 
localization of the fluctuations. If the phases are independent and uniformly distributed then the perturbations are 
delocalized. 

Finally there may be correlations between the various scales defined by i. In the language of |ljj this is what is 
called connectivity of the fluctuations. These correlations measure how much coherent interference is allowed between 
different scales, a phenomenon required for the rather abstract shapes and localization on each scale to become 
something visually recognizable as shapeful or localized. As in JjiJ this may be cast into inter-£ correlators. As we 
shall see these are in fact quite complicated for general anisotropic Gaussian theories. Therefore we have chosen not 
to dwell on this aspect of large-scale non-Gaussianity in this paper. 



III. A POSSIBLE METHOD FOR INTRODUCING GAUSSIAN FLUCTUATIONS IN ANISOTROPIC 

UNIVERSES 



We now present a po ssible way of introducing Gaussian fluctuations in anisotropic Universes such as the Bianchi 
models. In Section VII we will present another context in which anisotropy appears: periodic Universes. There we 
shall present more specific calculations of anisotropic Gaussian perturbations. Here we shall however use a method 
which relies simply on inspecting the reduced symmetry group anisotropic Gaussian perturbations must satisfy. This 
is a simple, if somewhat phenomenological, way of introducing the most general Gaussian perturbation which can live 
in an anisotropic background. Without actually performing a detailed perturbation analysis of these spacetimes, one 
can refine the analysis of || by using this prescription and possibly find more stringent constraints. 

Let an all-sky temperature anisotropy map be decomposed into spherical harmonics as in Eqn. (|lj). Then, for a 
general Gaussian theory, the a m are Gaussian random variables specified by a covariance matrix which must satisfy 
the symmetries of the underlying theory. In Friedman models the symmetry group is SO (3), but the symmetry group 
may be smaller. Anisotropic Gaussian fluctuations may be defined as Gaussian fluctuations with a covariance matrix 
satisfying a symmetry group which picks a favoured direction in the sky. We concentrate on anisotropic fluctuations 
with an SO (2) symmetry, that is, with cylindrical symmetry. 



2 It is curious to note that white noise has less structure than generic Gaussian fluctuations, but it also has more symmetry. 
It is tempting to associate reduction of symmetry and addition of structure. Anisotropic fluctuations have less symmetry than 
isotropic fluctuations, but they also have more structure, reflected in their non-Gaussian structure. 
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The general form of the covariance matrix may be obtained just from the representation theory of the symmetry 
group. The symmetry group breaks the {a l m } space into irreducible representations (Irreps). The a m may then be 
reexpressed in a basis adapted to these Irreps. Using Schur's Lemmas ||(| one knows (see Appendix I for more 
detail) that the covariance matrix of the theory must be a multiple of the identity within each Irrep []. Furthermore 
correlations between different a m can only occur for elements of different but equivalent Irreps. Hence, for any 
Gaussian theory subject to a symmetry which does not lead to equivalent Irreps, the spherical harmonic coefficients, 
expressed in a basis adapted to the partition into Irreps, must be independent random variables, and their variance 
must be a function only of the Irrep they belong to. As we shall see it may happen that the variance is the same for 
a set of Irreps. This degeneracy then leads to an accidental enlarged symmetry. If some of the Irreps are equivalent 
then in principle one may also have correlations between coefficients belonging to different but equivalent Irreps. 

As an example consider an isotropic theory. Then the {a^} for each £ are an Irrep of the symmetry group 5*0(3) 
represented by the D matrices 

iZMM)a^ = Dw(^0,#&, (12) 

where (?/>, 9, <p) are Euler angles. None of these Irreps is equivalent, as indeed none of them have the same dimension. 
Hence for a Gaussian isotropic theory the a l m must have a covariance matrix of the form 

( a m a m/) = Su'Smm'Ci (13) 

If the angular power spectrum Ce happens to be a constant (white-noise) over a certain section of the spectrum then 
this degeneracy increases the symmetry group of the theory: rotations among different £'s are now an extra symmetry. 
This is an accidental symmetry resulting from the degeneracy displayed by the particular model considered (white 
noise) and not required by the underlying theory. 

Now suppose that the symmetry group is SO(2), that is, the unperturbed model supporting the fluctuations is 
cylindrically symmetric. Then there is a favoured axis in the Universe and with respect to this axis the symmetry 
transformations are 

R{<p)a l m = e im <f-al (14) 

The Irreps are now indexed by £, m with m > 0. They are one dimensional complex Irreps for m > 0, and one 
dimensional real (and trivial) Irreps for m = 0. For the same m Irreps with different £ are equivalent Irreps. For each 
I we have a single Irrep of SO(3) which splits into I + 1 Irreps of SO(2). The covariance matrix of the theory now 
has the general form: 

{ a m a m}) = 8mm'Cfm\ (15) 

and we may call the diagonal terms Cg m of Cj^i the cylindrical power spectrum. It may now happen that = 
6 C| m |, and furthermore that a given model displays the degeneracy Cg\ m \ = Ci, that is the cylindrical power 
spectrum is white noise in m. In this case the 50(3) symmetry is accidentally restored. However this is no different 
from the white-noise model Ci — const referred to above. It is merely an accidental enlarged symmetry displayed by 
a concrete model and not a fundamental symmetry imposed by the underlying model. 

Accidental symmetries (eg. family symmetry in particle physics) are always regarded with horror. If they happen 
to exist, sooner or later a fundamental principle is sought which will promote them from accidental to fundamental 
symmetries. If they don't happen to exist a priori, such as in the case of fluctuations in anisotropic models, then 
better not postulate them in the first place. 



IV. NON-GAUSSIAN EFFECTS ON THE ANGULAR POWER SPECTRUM 

Gaussian anisotropic theories display many of the novelties present in non-Gaussian theories, such as the texture 
models considered in pl| , ^2[ . They trade their added predictivity in terms of non-Gaussian spectra for larger cosmic 



3 Schurs' Lemma only applies to finite dimensional representations, such as the ones offered by the a m . If one instead looks at 
the real space maps ST/T, then the representation space is S 2 . This is infinite dimensional, and indeed the covariance matrix 
of Gaussian theories is not diagonal, and is specified by the two-point correlation function C(9). 
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variance error bars in the angular power spectrum. Also the observed Cg may be correlated, a phenomenon called 
cosmic covariance and present in the texture models in (ll] jl^] . Cosmic covariance (or C e aliasing) induces great mess 
when comparing predicted and observed power spectra. Correlations allow for each observed power spectrum to have 
more structure than the average power spectrum. This may result in the average power spectrum corresponding to 
nothing that any observer ever sees. More subtle methods for predicting power spectra are then necessary. Two 
prescriptions are given in fl2|| . 

A. Cosmic variance surplus 

For a Gaussian isotropic theory the angular power spectrum 



has the variance 



Here we use the notation C £ to denote the random variable and Cg to denote its ensemble average. For a Gaussian 
anisotropic theory this variance is 

<7 2 {C £ ) = (21+1)2 zl C im ( 18 ) 
m——g 

If we define the average cylindrical power spectrum by 



Cg = — 

e 21 



1 - 

— E C 'm (19) 



171= - 



then 



> (20) 
v ' ~ 21 + 1 K ' 

It is a simple analysis exercise to prove this inequality and show that it is saturated only when Cg m = Cg, that is 
when the fluctuations are isotropic. 

Generally we may interpret this result as a reduction in the number of degrees of freedom in the x 2 induced by 
anisotropy. Suppose, for instance, that a theory is strongly anisotropic so that only a few m modes among the available 
21 + 1 contribute to the power spectrum Cg, for a given I. Then, effectively, the observed power spectrum C is the 
result of these few modes. Since these are still Gaussian variables the observed power spectrum is a \ 2 , but with an 
effective number of degrees of freedom equal to the number of predominant modes. If for example all the power is 
concentrate on the m — mode, then the C e is a Xx ■ If au the power is in a m > mode, the C l is a xi- 

We may use the ratio between the actual cosmic variance of the theory and its Gaussian prediction to quantify how 
anisotropic the fluctuations are. Quantitatively let us call anisotropy in the multipole I to the quantity 

^=#0=^y](%V (2D 



m— - 



which varies between Ai = 1 for isotropic theories to Ag = 21 + 1 for cylindrically symmetric multipoles (for which 
all the power is in the m = mode). 



G 



B. Cosmic covariance 



There are also correlations between different C l . For I ^ £' we have that 

cov(C l ,C e ') = 1 cov (l«™| 2 ' l«l'| 2 ) ( 22 ) 

m,m' 

For two (possibly correlated) complex Gaussian random variables zi and Z2 with uncorrelated real and imaginary 
parts, it can be shown that cov[|zi| 2 , \z2\ 2 ] = (ziz%) + (Z1Z2) , and so 

^■ PW ^ (23) 

where m in the summation runs from — mm.[£, £') to mm(£,£'). The off-diagonal elements (in £,£') in C„ therefore 
induce correlations among the various observed C e . A possible, but model dependent, way to do away with these 
correlations is to rotate the C l among themselves so as to diagonalize the covariance matrix (pf). These rotated C ( 
will then be independent, and so their average value is a good prediction for what each observer will see. Also, as 
shown in |p^| , in the rotated basis the cosmic variance error bars tend to be smaller and approach their Gaussian 
minimum. Therefore cosmic covariance, and larger cosmic variance error bars can be dealt with by means of this 
trick. However this trick does depend on each particular model, and is not a Universal prescription applicable to 
every model. 



V. THE NON-GAUSSIAN STRUCTURES EXHIBITED BY ANISOTROPIC GAUSSIAN THEORIES 

Anisotropic Gaussian theories also display non-Gaussian structure in the senses given at the end of Sec.||, that is 
they produce non-trivial non-Gaussian spectra. Here we shall find the most general type of isotropic non-Gaussian 
structure which can be mapped from these theories. 

We shall consider the anisotropic covariance matrix in more detail. Let the matrix C„ be split into its diagonal 
and its off-diagonal X~^ parts 

C5=*"'Q H +X5 (24) 
Then Cg\ m \, and so the bilinear form in the exponent of the Gaussian distribution 

FlO x.-ni, I -EE^^m I 125) 

m Li' 



M% = CZ'- 1 = - £fr (26) 

W|m| <-^W' 

and so the distribution factorizes into a factor which reveals the structure inside each multipole, and a factor which 
reveals correlations between different multipoles. We shall analyze these two factors in turn. 

Let's first assume that = 0. Repeating the transformation presented in Section || but using a covariance 

matrix of the form (|l^) one ends up with a rather complex distribution which has the form: 



ri( ,, -B: i .s[ r r l \ ii) = FU-''.U' nl )x\x^-y t (27) 



Unless C( m — Cg, the are not uniformly distributed. Also the C f will in general not be a xfe+i> ana - ^ ne function 
F(C e , B^) will not factorize. This means that not only will correlations exist between the B l m but the B^ will also 
be correlated with the angular power spectrum. The phases (j/ m on the other hand will still be uniformly distributed 
and independent. The phases tell us nothing about Gaussian anisotropic fluctuations. 

Hence anisotropic Gaussian fluctuations, when seen from the point of view of an isotropic formalism, are an example 
of delocalized shapeful fluctuations (explored in some detail in p9l ) . In the next two sections we will explore in more 
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detail the particular type of non-Gaussian effects which Gaussian anisotropic fluctuations may induce. The shapes 
exhibited by these theories are not the most general shapes, because there must be a scale transformation in the (r m 
which would render the B l m uniformly distributed again. Clearly not all shapes have this property. 

On top of this if ^ the distribution F(a l m ) does not factorize into factors which only depend on one t. 
Correlations between the different I will then appear, which in the language of (l9) amount to the emergence of 
connected structures: different scales are allowed to interfere constructively. In this paper we will not explore this side 
of the problem in depth. Nevertheless we have identified the non-Gaussian structures into which anisotropic Gaussian 
fluctuations are mapped. These are the delocalized shapeful (and possibly connected) structures defined in p9[ , or 
rather, a subclass thereof. 

We should note that although the {C e , B^ n , <j>f n } decomposition is not 50(3) invariant, the {C 1 ,!!^} already are 
SO(2) invariant |. Since the phases contain no information whatsoever on Gaussian anisotropic fluctuations they do 
not count as a device for making predictions in these theories (as much as one does not compute B l m for Gaussian 
isotropic theories). Hence the set of variables {C e , B^} is suitable for representing invariantly the most general form 
of non-Gaussian fluctuation which can be mapped from Gaussian anisotropic fluctuations. 



VI. GLOBALLY ANISOTROPIC UNIVERSES 



A useful set of models in which to explore these concepts are the homogeneous, anisotropic cosmologies, also know 
as the Bianchi models [Q. One can describe Bianchi cosmologies in terms of the metric 

gp.v = -n^ny + o 2 [exp(2/3)] A se^ef , (28) 

where n a is the normal to spatial hypersurfaces of homogeneity, a is the conformal scale factor, (3ab is a 3 matrix 
only dependent on cosmic time, t, and are invariant covector fields on the surfaces of homogeneity, which obey the 
commutation relations 

4„-e4 = Cr£ c e£e?. (29) 

The structure constants Cg C can be used to classify the different models. We shall focus on open or flat models 
which are asymptotically Friedman. These can be obtained by taking different limits of the type Vllh model which 
has structure constants 

C 2 m = = 1, C 2 21 = C 3 31 = Vh (30) 

It is convenient to define the parameter x — \Jhj{\ — f2o)j which determines the scale on which the principal axes of 
shear and rotation change orientation. By taking combinations of limits of f2 and x one can obtain Bianchi I, V and 
VII cosmologies. 

We are interested in large-scale anisotropies so it suffices to evaluate the peculiar redshift a photon will feel from 
the epoch of last scattering (Is) until now (0) 

AT A (r) = (f l Ul ) - {f l Ul ) ls - / Vf k a jk dT (31) 

J Is 

where f = (cos sin 0, sin 9 sin 0, cost/)) is the direction vector of the incoming null geodesic, u is the spatial part of 
the fluid four- velocity vector and to first order, the shear is cry = d T (3ij. To evaluate expression (|3ll), one must first 
of all determine a parameterization of geodesies on this spacetime. This is given by 

tan(^i) = tan(y ) exp[-(r - T )Vh] 
0(t) = 8 + (t- t ) 

~ ln{sin 2 (^) + cos 2 (^) exp[2(r - r Q )Vh}} (32) 



4 We are assuming that not only the Universe is anisotropic but that we know, a priori, what its symmetry axis is, eg: by the 
detection of a Hubble-size coherent magnetic field. Alternatively we leave the Euler angles of this axis free, to be estimated by 
some MLE. 
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Solving Einstein's equations (and assuming that matter is a pressureless fluid) one can determine u and cry . A general 
expression for ( |3l| ) was determined in S: 

AT A (r) = (-^r) 2 ^ /1 o -^ x {[sin<focos0 o - sin <p ls cos 9i s (l + zi s )] 
Id s Zq 

T ° 3M V" 0) sin20[cos(0) + Bin(g)] , f T } (33) 
is "0 smh (vhT/2) 

A useful discussion of the different CMB patterns imprinted by the unperturbed anisotropic expansion is presented 
m j|. The patterns can be roughly said to be constructed out of two ingredients: a focusing of the quadrupole when 
fZ < 1 and a spiral pattern when x is finite. The Bianchi VII/j is most general form of homogeneous, anisotropic 
universes in an fZ < 1 which are asymptotically Friedman-Robertson- Walker. The pattern is of the form: 

^L = f x (e)cas{<i>-m) (34) 

In each 9 = const circle the pattern has a dependence in (f> of the form cos(0 — (j>). The phase (j> depends on 9, and 
hence the spiralling of the simple cold and hot bump induced by the cos0 dependence. The functions fi(9) and 
4>{9) are rather complicated functions which have to be evaluated numerically, and depend on various details of the 
particular Bianchi model within the type we have chosen. It is curious to note, however, that only the power spectrum 
Ct and the phases <f> are sensitive to these details. All the spirals imprinted by Bianchi VII^ models have moduli of 
the form 



pi l = S ml f 2 (l,x) (35) 



Therefore their shape spectra will always be 



B m = 1 for 2<m<£ 

B l m = for m = 1 (36) 

The background patterns in Bianchi VII^ models are all localized, shapeful, and connected structures. Depending 
on the model they will however have different positions, power spectra, and connectivity. Nevertheless, their shape 
spectra is always the same exact shape, of form (|3^), without any cosmic variance error bars. Confusion with a 
Gaussian is zero. Confusion with the shape of a perfect texture hot spot is zero as well. These have a non-Gaussian 
spectrum of the form 

B e m — 1 for 1 < m < £ 

(37) 

Although the shape spectrum is the same up to the last B^, the confusion between the two theories is zero. 

VII. AN ANISOTROPIC GRISCHUK-ZEL'DOVICH EFFECT 

We shall now consider an example of a flat homogeneous and isotropic universe with topological identification along 
one axis. This example is simpler than most of those considered in the literature but illustrates one of the key features 
of such models: the breaking of statistical isotropy in the fluctuations. Let us consider a universe with a topological 
identification along the z axis. All functions defined on such a space satisfy: 

$(x,y,z) = <P{x,y,z + L) (38) 

By considering a flat universe we can restrict ourselves to calculating the Sachs- Wolfe effect. The temperature 
anisotropy from the surface of last scattering will be given be: 

f- (n)= -^/ d2fc £ (39) 

j=-oa 

where A = r)o — rji s is radius of the surface of last scatter and q = (k cos </>, fcsin^, 27r-jj ) . We can expand the exponential 
in spherical harmonics to get 
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The a( m s are given by 

ae m = J fk^^Y^M (40) 

We now assume statistical homogeneity and isotropy of 8 and a scale invariant power spectrum 



(Stjfaj,) = 8 2 (k- k'^q- 1 q = X /k2 + (41) 



This leads to the covariance matrix for the a£ m s 



d ^3 Y em(q) Y t'm'W ( 42 ) 



Expressing Yi rn (n) — Pi m (cos9)e" n ^ and performing the azimuthal integral, one immediately finds that the covari- 
ance matrix is diagonal in m, so one has 



{a} m a PJm 



— —5 mm > kdk ^ 5 Ptm(q)Pi'm> (<?) (43) 

3 

It is convenient to define x — Afc an d Mi — "2iTaj where a = A/L, i.e. the ration of our horizon to the topological 



identification scale. If we now define y — Mj/y M 2 + X 2 the expression simplifies to 



(oJ m o^ m '> = -2 V 5mm ' / d 2/Z^ — —Ptm{y)Pt'm>{y) (44) 

4C j ^ 

In the limit where the identification scale goes to infinity we get the standard result 

Z" 1 1 1 

lim {a} m a£' m >) oc / dyje(-)5 u >5 mm > cx -^8w5 mm ' (45) 

i.e. a scale invariant, diagonal covariance matrix. In the case of finite a this is not the case. Consider the quadrupole. 
The ring spectra has two components, B\ and B 2 with a probability distribution function 

F(r ' Bl ' B2) = C Z%'^ X ex P^ 2 /2<x 2 2 )[c 2 3 2 2 (l + ^B 2/3 )]} (46) 
where we have defined 

o\ = «|a ao | 2 )<|o ai | 2 )<|a a2 | a )) 1 /3 

By exploring the dependence of ci and c 2 on a we can see how the probability distribution function of the Bs change 
with topology scale; to a very good approximation we find 

Ci = y/2 

c 2 2 

2a .,„, 
c 2 = — — 48 
1 + a 

We can see the signature for non-Gaussianity arising here. For a non zero a there are correlations between the three 
statistical quantities, in particular (B1B2) cx a/(l + a). The probability distribution function for B\ and £? 2 (defined 
on [0, 1] x [0, 1]) becomes peaked at 0. We have focused on the quadrupole where the effect is easy to see. The method 
is systematic however, and one can construct the probability distribution function of the high order ring spectra in 
the same way. 

This a curious application of the idea put forward in plj ], an anisotropic Grischuk-Zel'Dovich effect. By looking at 
the shape of the low t multipolc moments we can constrain the degree of statistical anisotropy outside the current 
horizon. Note that, already for a < 1 there are deformations in the covariance matrix which may be statistically 
significant with current data. This will be pursued in a future publication p2l . 
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VIII. CONCLUSIONS 



In this paper we have presented a new technique for quantifying non-Gaussianity on large scales. It is the extension 
of the non-Gaussian spectra developed in jl9| to the surface of the two sphere. As we have shown in Section [n] the 
construction is slightly different to take into account the particularities of the spherical harmonic basis. However 
the quantitative interpretation of the different levels of non-Gaussianity follows through, exactly as in |l9| . One can 
identify the information contained in the ring, interring and phase spectra with shape, connectivity and localization. 

An interesting and untapped application is to universes with statistically anisotropic fluctuations. Developing the 
idea put forward in (hJ we explain how statistical anisotropy and non-Gaussianity are intimately related. From this 
one can infer some novel properties of the covariance matrix of fluctuations in statistically anisotropic spacetimes. 
In particular, features which appear in non-gaussian theories of structure formation, like textures [pj],I2| will appear 
here: a surplus of cosmic variance and cosmic co-variance of the power spectra. 

There are a number of situations where these results are applicable. One is in the case of anisotropic universes, i.e. 
universes which aren't Friedman-Robertson- Walker. There are a number of known examples 0,^]- Without actually 
doing perturbation theory on them we argue for a natural prescription for adding fluctuations in the CMB to such 
models. It consists of finding the reduced symmetry group of the temperature patterns and constructing anisotropic 
Gaussian random fields with such properties (in which the covariance matrix satisfies those symmetries, and not 
more). In fact it can be shown that these symmetries can be deduced from the geodesic structure of the space time 
p3| . This can be a first step in extending the prescription used in |(| for constraining general anisotropic models with 
the COBE 4 year data. A brief analysis is made of the relevant Bianchi models for which we present the non-Gaussian 
spectra. 

Another, different application is the case of homogeneous isotropic models where an anisotropic topological iden- 
tification has been imposed. As an example we identify one direction in space. One finds that statistical isotropy 
is broken. This can be easily from the following: if we look along the axis of identification, and the identification 
scale is smaller than our horizon, one will find strong correlations between patches of the microwave sky which are 
reflected about the uncompactified plane 0. By looking at the structure of the covariance matrix one can see that 
this anisotropy will manifest itself by inducing not only non-Gaussian ring spectra but also inter-ring spectra. This 
non-Gaussian manifestation may persist if we consider the identification scale to be large than our horizon. We name 
this effect the anisotropic Grischuk-Zel'Dovich effect. 

Although we now have a high quality measurement of anisotropies on large angular scales we are confronted with 
the hardships of the real world. Galactic contamination leads one to consider an anisotropic rendition of the sky and 
considerably complicates the analysis of the COBE four year data. It is well known that one of the consequences is that 
the quadrupolc measurement should be viewed with scepticism. Unfortunately it is the quadrupole which could supply 
us with a probe of primordial CMB on the largest angular scales. There may be ways around these shortcomings. 
One can try and reformulate our non-Gaussian spectra on the largest angular scales using the techniques put forward 
in p4| . This would involve a proper likelihood analysis and to make the problem tractable one would need to find 
an adequate parametrization of the non-Gaussian spectra. The fact that we have devised a consistent method for 
characterizing non-Gaussianity on all scales may allow us to use the cumulative information of all I > 2 to infer 
the behaviour on large scales; all modes will be affected to some extent by large scale anisotropy. Finally it would 
be interesting to analyse in more detail the observability of the anisotropic Grischuk-Zel'Dovich effect, taking into 
consideration issues of cosmic variance. 
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APPENDIX I - GAUSSIAN FLUCTUATIONS WITH A SYMMETRY 



In this Appendix we detail the group theory argument sketched in Section III. This argument is not necessary in 



the SO(2) case targeted in this paper, where the covariance matrix may be easily derived directly. However, it opens 
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doors to more general anisotropic symmetry groups. It also shows how the general form of the covariance matrix 
depends not on the symmetry group, but only on its irreducible representations (Irreps). 

Fluctuations (Gaussian or not) in any Universe must be subject to the symmetries of the underlying cosmological 
model. However, due to the random nature of the fluctuations, they must satisfy these symmetries only statistically. 
By this we mean that the statistical ensemble of fluctuations, and not each realization, should be subject to the 
symmetries. If a symmetry transformation is applied to each member of the ensemble, then each member may 
change, but the ensemble should remain the same. For instance, for the 50(3) symmetry group only realizations 
containing only a monopole £ = are left unchanged by rotations. Nevertheless much more general fluctuations respect 
statistical isotropy. Similarly, only the m = modes are cylindrically symmetric, but more general fluctuations are 
50(2) statistically invariant^]. 

Gaussian fluctuations are fully specified by their covariance matrix 

^ / m 1 m 2 ~ rmi"m2/ 

which may be seen as a bilinear form on the {a^} space. Hence the statistical symmetries of a Gaussian theory are 
equivalent to the requirement that the covariance matrix is left unchanged by any symmetry transformation. Let 
G be the symmetry group of the underlying cosmological model as projected on the sky. Let's first suppose that G 
breaks the {a^} into a set on non-equivalent Irreps. Then, let's find a new basis a M adapted to G, where L now 
labels the Irrep the basis element belongs to, and M the actual element. Then G is represented by a set of matrices 
G\, [M , acting on a\ { as 



a-M — Ga M — G\ AM ,a M , (50) 



Cm\m 2 = «4?*> (51) 



The covariance matrix for the a M 



must remain unchanged by the transformation (|50|), so that 

< ^M 1 M 2 — \ a M 1 a M 2 I — Mi M{ M 2 M[ M' 2 ~ Ly M 1 M 2 \ oz ) 

which for unitary representations amounts to the commutation relation 

G Ll C LlL2 = C LlL2 G L2 (53) 

Let us now recall Schur's Lemmas pOj . 

Schur's Lemma 1 Let T and V be two Irreps of a group G with dimensions d and d! , and let there be a dx d' matrix 
A such that 

T(g)A = AT'(g) (54) 
for all group elements g. Then either A = or d = d' and detA ^ 0. 
It follows that if A ^ then T and T' are equivalent. 

Schur's Lemma 2IfT is a d-dimensional Irrep of a group G and B is a d x d matrix such that 

T(g)B = BT(g) (55) 

for all group elements g, then B = XI. 



5 Some texture models exhibit approximate SO(2) symmetry at low £ in each realization. Then an axes system exists in which 
the m = mode has much more power than any other. This is of course a very non-Gaussian effect which cannot be simply 
reproduced by anisotropic Gaussian fluctuations. 
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Combining the two Lemmas we can then find that the covariance matrix must take the form 

Cm\m 2 = (vmSmD = S LlL2 S MlM2 C Ll (56) 

An interesting result is that if G breaks the {a^} space into non-equivalent Irreps spanned by {aj^}, then these must 
be independent random variables with a variance which can only depend on the Irrep they belong to. This argument 
applies for instance if the symmetry group is SO(3), in which case L = £ and M = m. 

The argument just present breaks down however, if some of the Irreps defined by G are equivalent. Let the {&„} 
space now be spanned by an adapted basis {aj^f}, where L labels each class of equivalent representations, D labels 
the actual Irrep within this class, and M labels the elements within each Irrep. Then, from the second Schur's Lemma 
we know that within the same Irrep we still have: 

{a L M D 1 a L M D 2 *)=5 Ml M 2 C LD (57) 

but for different but equivalent Irreps (D\ ^ D2) we now have 

(<F:< 2 *)=C L AI X 2 (58) 
with detC 7^ 0, whereas for non-equivalent Irreps {L\ ^ L2) we still have 

(*M 1 Dla m D2 *)=° ( 59 ) 

Therefore, although the are independent within each Irrep and among different non-equivalent Irreps, correlations 
may exist between different but equivalent Irreps. Of course one may always rotate the a\P within each class of 
equivalent Irreps so as to diagonalize the covariance matrix. However such a rotation is model dependent, and cannot 
be determined from the symmetries. This situation happens for instance in the case of SO(2), with L — m and D = £ 
(no M index, since the Irreps are one dimensional). Each m provides a class of equivalent representations, with the 
same m but different £. Then the covariance matrix takes the general form 

("Vm*) = &m 1 m 2 C t t ^ 2 (60) 

and as we see although correlations among different m are not allowed, now we may have correlations between different 
£, for the same m. We could rotate the a m in I for each fixed in so as to diagonalize the covariance matrix, but such 
procedure naturally would depend on the covariance matrix one starts from, and would therefore be model dependent. 
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